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Abstract 

A step towards a successful implementation of time- 
resolved digital in-line holography with extreme ul- 
traviolet radiation is presented. Ultrashort XUV 
pulses are produced as high-order harmonics of a fem- 
tosecond laser and a Schwarzschild objective is used 
to focus harmonic radiation at 38 nm and to produce 
a strongly divergent reference beam for holographic 
recording. Experimental holograms of thin wires are 
recorded and the objects reconstructed. Descriptions 
of the simulation and reconstruction theory and algo- 
rithms are also given. Spatial resolution of few hun- 
dreds of nm is potentially achievable, and micrometer 
resolution range is demonstrated. 
PACS 42.65.Ky; 42.40.-i 

1 Introduction 

Advanced coherent and short-pulse photon sources in 
the XUV and X-ray range are being extensively de- 
veloped with the goal to provide both temporal and 
spatial resolution, thus allowing scientists to explore 
at the same time structure and dynamics of matter. 
Examples of this development are the effort to build 
X-ray free electron lasers at various places in the 
world [1] and to develop seeded soft-x-ray lasers [2]. 
One of the first experiments performed at FLASH, 
XUV-free electron laser facility in Hamburg is time- 
resolved X-ray holography using a clever single-shot 
design [3]. Much smaller sources are those provided 
by high-order harmonic generation in gases [4] . These 
sources have the advantage of being table-top and 
affordable in university laboratories. The temporal 
resolution reaches the attosecond time scale, while 
the wavelength range varies from the UV to the soft 



X-ray spectral regions. The highest conversion effi- 
ciencies, of the order of a few times 10 -5 , and energy 
per pulse (fractions of /iJ) are obtained when Xe or 
Ar are used for generating the harmonics, with wave- 
lengths typically between 30 and 50 nm. The spatial 
coherence has been measured to be excellent, compa- 
rable to that of the generating laser beam [5,6]. 

These properties make this radiation an interesting 
candidate for time-resolved microscopy [7], in par- 
ticular using holographic approaches [8], at least in 
the long- wavelength part of the XUV range (around 
30 nm) where the number of photons per pulse and 
per harmonic can be as high as 10 9 — 10 10 . A few 
promising experiments have been performed, using, 
in particular, digital holographic techniques in the 
XUV [5,9] and X-ray range [10-13]. Bartels et al. [5] 
used a very simple setup, where an object was placed 
directly in a harmonic beam, containing harmonics 
of orders 17 through 23. The resulting interference 
patterns, produced by the diverging harmonic beam, 
acting as a reference wave, and the waves diffracted 
from the object, were recorded by an X-ray CCD 
camera. The resolution achieved was not very high 
due to a low numerical aperture, and the use of sev- 
eral harmonics blurred the image. However, images 
of the tip of a scanning microscope could be recon- 
structed. Morlens et al. [9] used two multilayer spher- 
ical mirrors to focus the beam on the sample in or- 
der to achieve a larger numerical aperture and to se- 
lect only one harmonic, the 25 th . Holograms of spi- 
der threads were recorded and reconstructions with 
resolution below one micrometer could be success- 
fully performed. In both cases, however, acquisitions 
of over 1000 shots were needed to achieve high sig- 
nal/noise ratio. Recently Tobey et al. [13] performed 
time-resolved holography using high-order harmon- 



ics. The experiment was done with a simple single- 
reflection geometry and without any focusing optics. 
In the present work, we use high-order harmonics pro- 
duced by a powerful 10 Hz, 35 fs, titanium-sapphire 
laser system to perform digital in-line holography of 
micrometer wires at 38 nm wavelength, using focus- 
ing optics with higher spatial resolution than in the 
previous experiments. High quality holograms are 
obtained by adding about 100 images. Single shot 
images are also obtained, though with low signal to 
noise ratio. Images are successfully reconstructed by 
using a reconstruction algorithm. In Section 2 we de- 
scribe the experimental method and results, while in 
Section 3, we give details of our reconstruction algo- 
rithms and apply them to the analysis of the exper- 
imental holograms. In Section 4, we discuss how to 
improve the experiments to realise high-quality sin- 
gle shot digital in-line holography, leading the way 
to time-resolved holographic microscopy in the fem- 
tosecond regime. 

2 Experiment 

Fig. [I] describes our experimental setup. The laser 
beam is a low-energy arm of the 40 TW laser sys- 
tem at the Lund Laser Center, a titanium-sapphire 
laser based on the chirped-pulse amplification tech- 
nique. Pulses are delivered at a 10 Hz repetition rate, 
35 fs duration, 800 nm central wavelength and up to 
100 mJ energy. 10 mJ are actually used to gener- 
ate harmonics by focusing loosely (1 m focal length) 
the laser beam in an argon gas jet. After that the 
fundamental beam is blocked by a 200 nm-thick alu- 
minium filter. Fig. [2] shows a typical harmonic spec- 
trum recorded with a flat-field XUV spectrometer. 
A multi-layer coated Schwarzschild objective is used 
to focus tightly the harmonic beam and at the same 
time to select a narrow spectral region. An object 
(wire) is located a few mm after the image point of 
the harmonic source and holograms are recorded by 
using imaging multichannel plates (MCP) coupled to 
a phosphor screen. A CCD camera is used to record 
the image on the phosphor screen. 

A Schwarzschild objective is an optical system con- 
sisting of two spherical mirrors, one convex and one 
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Figure 1: Experimental setup. The laser beam is focused into 
a gas jet producing harmonics, which then pass trough an alu- 
minium filter. A multi-layer coated Schwarzschild objective 
focus the harmonic beam on an object and the resulting in- 
terference pattern is recorded by a multi-channel plate and a 
CCD camera. 



concave of which the centers of curvature coincide. 
A tight focusing can be obtained. In our design the 
two radii of curvature are equal to 3.39 and 9.17 cm 
respectively, with a resulting focal length of 2.69 cm. 
The Schwarzschild objective is placed 200 cm from 
the gas jet producing the harmonics and the detec- 
tor 30 cm further away. The object is located about 
3 mm after the focus of the Schwarzschild objective. 
To avoid loosing significant XUV energy, we use it 
in a slightly off-axis geometry as shown in Fig. [I] 
without introducing significant aberration. Due to 
this off-axis configuration, the beam going out of the 
objective is forming an angle of 10° with the ingoing 
beam. It is coated with a multilayer structure giving 
a transmission window centered at A = 37 ± 3 nm. 
The maximum transmission is approximately 40 %. 

The wavelength transmission of the objective al- 
lows us to select the 21 st harmonic. The inset at 
the top of Fig. [2] represents the far-field profile of 
that harmonic, with a 1 mrad divergence. The en- 
ergy available in the 21 st harmonic is measured to be 
5 nJ which corresponds to 10 9 photons per pulse, re- 
duced by a factor of 5 after transmission through the 
Al filter. 

The theoretical resolution of the experiment is 
given by: 

R _ A _XL 

2tan(9) d ' K ' 

where 9 is the half angle of aperture, A the wave- 
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Figure 2: High harmonic spectrum and profile of the 21 st har- 
monic. The dashed line represents the transmission of the 
Schwarzschild objective. 




Figure 3: Holograms recorded after 200 shots with different 
objects: a hair of 65 /im diameter (a), a wire of 18 diameter 
(b) and a spider thread of 5 /im diameter (c). 



length of the radiation used, L the distance between 
the object and the detector and d the diameter of the 
detector area illuminated by the beam. In our setup, 
the half angle of aperture after the Schwarzschild ob- 
jective is 4.3°, leading to a resolution of 250 nm. 
Recorded holograms of wires of diameters varying 
from a few tens to a few micrometers are shown in 
Fig. [3J These results are obtained after integration 
of 200 images. 



3 Reconstruction 

In this Section the algorithm used for simulation of 
holograms and reconstruction of experimental images 
is presented. Many innovative techniques have been 
put together to allow reconstruction of experimental 
datas. Due to the tight focusing of the beam on the 
object, the beam has gaussian properties, that have 
been taken into account. Aliasing is also avoided in 
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Figure 4: Illustration of the theoretical formalism and nota- 
tions. Two Fresnel planes, the object plane (O) and the holo- 
gram plane (H), are indicated. 



this process. Furthermore the problem of the twin 
images is solved by introducing an iteration process. 

Our reconstruction algorithm is based on the 
Fresnel- Kir choff integral, illustrated in Fig. [4] [14]. 
Every point in the object plane O (x Ql y Q ) is the 
source of a spherical wavelet. Their addition in the 
hologram plane H (xh, y^) is the electromagnetic 
field : 



E H (x h: y h ,Z + L) = 



i f e ~ikr 

- / E (x ,y ,Z) 

A Js r 

x cos 0dx o dy O: (2) 



where r is the distance from the point Pi to the point 
P2, is the angle between P1P2 and the normal to the 
planes, S is the area over which the integral is carried 
out. Z and L define the position of the object and 
hologram planes as illustrated in Fig. [4j 

In the paraxial approximation, \xh~ x Q \, \yh~Vo\ <^ 
L, cos# r can be expressed to first order as: 



L- 
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(3) 



Using 



Fk{x 



y Q ,Z) = E (x ,y ,Z)exp(-ik^j^) 



and its Fourier transform Fq, the resulting field can 
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Figure 5: Far-field, object and hologram planes, F, O and H, 
respectively. A new hologram plane (H') closer to the object 
plane is also introduced to solve the problem of the pixel size 
(see text). 



be written as 



E H (x h ,y h ,Z 
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(4) 



We now consider a third plane, (F), called far-field 
plane, far from the object and hologram planes, as 
shown in Fig. [5] The electric field in the far-field 
plane can be determined from the electric field in any 
of the two other planes [15], as 
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L\ and are the distances between F and O, and 
between F and H, respectively (see Fig. [5^. By 
dividing the two equations an expression for F^ 2 as 
a function of Fq 1 is obtained, and vice versa. 
To simulate a hologram we can then calculate 
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Similarly, to reconstruct the object, we use 

E (x ,y ,Z) =e ik *-^r 



FT~ 1 [G'Ffl 2 ( 
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where 
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This reconstruction method works quite well for 
plane or collimated waves. However in our case, tight 
focusing optics is placed before the object, leading 
to a strong beam divergence. The beam cannot be 
approximated as a plane wave anymore. Moreover, 
the typical size of the object is orders of magnitude 
smaller than the size of the detector and the pixel 
sizes in the hologram and far-field planes are large. 
The phase factors of different terms of Eh and Eo 
therefore vary very rapidly, leading to aliasing and 
picture replicas [16]. To overcome this problem, we 
introduce a new hologram plane, hologram plane (H f ) 
(see Fig.|5|. This plane is situated close to the object 
plane and there is almost no magnification between 
the two planes so that the algorithm described above 
can be applied. A simple transformation is then made 
to obtain the hologram in H . More formally, the field 
in the object plane can be written as 

-\- 

Eo(x ,y ,Z) = E (x ,y ,Z)exp(-ik ° ° ), 

(10) 

where Eo is a slowly varying term (product of the 
object transmission field and the slowly varying part 
of the reference field, taken to be Gaussian) while the 
phase term is a rapidly varying spherical term. We 
apply the following coordinate reduction: 

U u 

(Xh,Vh) (Xh'iVh') = (-jrXh, -jrVh), (11) 

with L' = LZ/(L + Z). After some elementary 
mathematical manipulation, the field in the hologram 
plane can be expressed as: 



E H (x h ,y h ,Z + L) = 



Z-V 



-ik(L-L') 



-ik(x 2 h ,+y 2 h ,) 



f) ^E H ,(x h ,,y h ,,Z + L f ).(12) 
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where En'ixh' ,yh' ? Z + V) is the field in H' corre- 
sponding to the slowly-varying field Eq in the ob- 
ject plane. The terms in front Eh> in Eq. (12) are, 



apart from a constant factor, pure phase terms that 
vanish when calculating the hologram intensity. This 
equation can be interpreted as follows: the diffraction 
calculation with a very divergent reference wave is 
equivalent to the diffraction calculation with a plane 
wave, on a plane H' separated from the object by 
L', followed by a magnification by (L + Z)/Z. The 
simulation algorithms described before can now be 
applied on Eq, using H' . To get the real hologram 



in H we just make a scale change using Eq. ( 11 ). The 
reconstruction algorithm is based on the same idea. 

In the reconstruction process, the phase informa- 
tion is lost because only the intensity of the holo- 
gram is recorded. To retrieve this information, an it- 
eration process inspired from the Gerchberg-Saxton 
algorithm, similar to the one described in [17], has 
been implemented. The aim of this process is to 
eliminate the twin images following reconstruction in 
in-line holography. Here are the different steps: 

1. Reconstruction of the object, using the intensity 
information of the hologram (no phase informa- 
tion). 

2. Simulation of a new hologram using the real part 
of the reconstructed object as object field. 

3. Replacement of the intensity of this new holo- 
gram by the intensity of the original hologram. 
The phase of the simulated hologram is now used 
as a first guess of the phase of the field on the 
hologram plane. 

4. Reconstruction of the object using this last field 
guess. 

5. Repeat steps 2 through 4 until convergence. 

Fig. [6] illustrates the performance of this algorithm 
on simulated data. In (a) the object is represented 
and in (b) the simulated hologram. The object is re- 
constructed in (c) without iteration and in (d) with 
100 iterations. The iteration process clearly improves 
the quality of the reconstruction. We also studied the 
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Figure 6: A hologram is simulated in (b) from the object in (a) 
with one harmonic. The object is reconstructed in without (c) 
and with iterations (d). In (e) the hologram is simulated with 
three harmonics and reconstructed in (f). The holograms are 
simulated and reconstructed with the experimental conditions: 
L = 0.27 m; Z = 1 mm; pixhologram — 95 /im; L\ = 10 6 m. 
The object is illuminated by a Gaussian beam with a half angle 
of aperture of 4.3 °. 



simulation and reconstruction of an hologram gener- 
ated by many harmonics. In (e) the hologram is sim- 
ulated with three harmonics (the 19 th , the 21 st and 
the 23 rd ). Despite the relative blurring of the fringes 
in the hologram, the object can still be reconstructed 
in (f) also with 100 iterations. 

In the experiment, the images are recorded by a 
digital camera and processed by a computer program. 
Thus the influence of the hologram pixel size has to 
be taken into account to estimate the experimental 
resolution, which will depend on the pixel size of the 
reconstructed object given by the algorithm. As seen 
in Eq. ([6| and ([sj) we need to calculate numerically 
several 2D Fourier transforms and their inverses to 
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simulate holograms or reconstruct object from holo- 
grams. These calculations involve a discretization of 
H and O in arrays of pixels. Understanding the rela- 
tions between these pixel sizes is crucial. The size of 
the reconstructed pixel of O has indeed to be of the 
order of the theoretical optical resolution in order to 
reach the full power of the method. Numerically, a 
reconstructed object with pixel size smaller than the 
optical resolution will exhibit diffraction blurring of 
the order of the optical resolution. For simplicity we 
consider the shape of the pixels of the object and of 
the hologram as squares. 

In the case of Fresnel diffraction calculated with a 
2D Fast Fourier Transform, the pixel sizes in F and 
in O are related by the simple expression: 



PIX far field 



XLi 



1 



(13) 



N pix object 

where L\ is the distance between the far field plane 
and the object, N is the number of pixels in one raw 
of the object plane, pixf ar fi e id is the length of the 
pixel side in F and pix bj ect is the length of the pixel 
side in O. A similar expression relates the pixel size 
in H and in F, so that pix ho i ogram = pix object L 2 / / L 1 . 
Taking into account the reduction factor introduced 
before, we have 



object 



L X Z 



Li(L- 

Z 
L + Z 



Z) - LZ 



P^X hologram 



piXhol 



ogram • 



(14) 



The pixel size of the reconstructed object depends on 
Li, L and Z and on the pixel size of the hologram. 
The minimum value is reached when L\ — > oo. To ob- 
tain an even smaller pixel size, we should thus reduce 
the pixel size of the hologram by interpolating the ex- 
perimental hologram or recording with a device with 
smaller pixels. Note however that, the resolution of 
the experiment is not limited by the reconstruction 
algorithm but only by the optical resolution given by 
Eq. Q. 

The recorded holograms are reconstructed using 
the algorithm described above and are presented in 
Fig. [7| We observe that the reconstruction is pos- 
sible for all the objects and that the phase retrieval 




Figure 7: Holograms recorded after 200 shots with different 
objects: a hair (a), a wire of 18 jum diameter (d) and a spider 
thread (g). Reconstruction of the holograms without iteration 
(b, e, h) and after 100 iterations (c, f, i). 




Figure 8: Hologram recorded in a single shot of a wire of 18 
[iu\ diameter in (a). Reconstruction of the hologram without 
iteration (b) and after 100 iterations (c). 



process is efficient. It is possible to record and recon- 
struct holograms of wires of 6 /im diameter, which 
is still far above the theoretical resolution of our ex- 
perimental setup (250 nm). By comparison a single 
shot image of a 18 jam wide wire is shown in Fig. [8j 
The reconstruction is still possible but with reduced 
quality, due to the low signal-to-noise ratio. 



4 Conclusion 

Holograms of thin wires are recorded using high or- 
der harmonics at 38 nm. A good spatial resolution is 
achieved with the help of a Schwarzschild objective. 
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Since single shot holograms are recorded and recon- 
structed, we are able to reach the femtosecond tem- 
poral resolution range. A reconstruction program for 
digital in-line holography has been developed, tested 
and implemented. An iteration process was also in- 
troduced to improve the phase retrieval process of the 
reconstruction. With the help of this program we are 
able to reconstruct different objects. 

This work proves that digital in-line holography 
using XUV radiation is a possible alternative to mi- 
croscopy and a step in its optimization has been 
taken. Although there are many improvement that 
can - and should - be done, this is a step towards suc- 
cessfully performing time-resolved holography with 
extreme ultraviolet radiation. The main improve- 
ments that will be made are mostly experimental. 
It seems important to select only one harmonic or- 
der to have a hologram which can be entirely recon- 
structed. This could be possible by adding other mul- 
tilayer optics efficiently selecting one harmonic order. 
Simultaneously, the harmonic photon flux should be 
improved by at least a factor of 10 to increase the 
signal-to-noise ratio. Promising techniques using, for 
example, improved phase matching [18, 19] or multi- 
color fields [20,21] make this within reach. In the fu- 
ture the possibility of single-shot measurements gives 
this technique a great potential for imaging femtosec- 
ond processes. 
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